EPJ manuscript No. 

(will be inserted by the editor) 



Similar final states from different initial states using 
new exact solutions of relativistic hydrodynamics 

M. Csanad^'^, M. I. Nagy^^'^, and T. Csorgo^-^ 

^ ELTE, Eotvos Lorand University, H - 1117 Budapest, Pazmany P. s. 1/A, Hungary 
2 MTA KFKI RMKI, H-1525 Budapest 114, FOBox 49, Hungary 

Abstract. We present exact, analytic and simple solutions of relativistic perfect 
fluid hydrodynamics. The solutions allow us to calculate the rapidity distribution 
of the particles produced at the freeze-out, and fit them to the measured rapidity 
distribution data. We also give an advanced estimation of the energy density 
reached in heavy ion collisions, and an improved estimation of the life-time of the 
reaction. 

1 Introduction 

Since the birth of hydrodynamics, there were many interesting exact solutions found to the 
complicated non-linear coupled differential equations of them. In contemporary research, such 
as the description of collective properties of high-energy elementary particle and heavy ion 
reactions, one has to deal with not only hydrodynamics but relativistic hydrodynamics. This 
is an even more complicated topic, and likely this is the reason why only a few exact solutions 
exist for relativistic hydrodynamics, in contrast to the nonrelativistic case. 

Exact solutions are important in at least two ways: first, they can be used to test numerical 
codes reliably. Second, they provide an invaluable insight into the details of the evolution of 
the matter created in high-energy reactions. After recalling some presently known solutions, 
such as the Landau-Khalatnikov solution and the Hwa-Bjorken solution, we present a class 
of accelerating exact and explicit solutions of relativistic hydrodynamics. These solutions are 
advantageous compared to the previously mentioned ones. Then we show how the new solutions 
can be applied to describe the evolution of the matter created in high energy reactions. Our 
treatment uses simple formulas and takes the presence of acceleration into account. 

This article is based on refs. |1|2|3| . These results are extended here with more detailed 
simulations that indicate that even for a fixed equation of state, different initial conditions 
could lead to the same freeze-out distributions at mid-rapidity. However, the width of the 
rapidity distribution can be utilized to select from among these time evolution scenarios. 

2 Basic equations 

In this section we recapitulate the equations of relativistic hydrodynamics. We will use the 
following notations: = (t, r) is the coordinate four-vector, r = {rx,ry,rz) is the coordinate 
three-vector. The metric tensor is g^j^ = diag {1, —1, —1, —1), the four-velocity field is u'^, 
the normalization is u'^u^ = 1. The three-velocity v is defined as u'^ — 7(l,v), with 7 = 

e-mail: csanad@elte.hu 
^ e-mail: nagymarciSrmki . kf ki . hu 
e-mail: csorgo@sunserv.kfki.hu 



2 



Will be inserted by the editor 



(l — v^) ^ ^ . We denote the so-called conioving derivative by that is, ^ = ^ + vV. The 
thermodynamical quantities are the following: s is the energy density, p is the pressure, w — e+p 
is the enthalpy density, T is the temperature, a is the entropy density. When there are some 
(conserved or non-conserved) charges present, we denote them by n^, and the corresponding 
chemical potentials by /i^. 

The equations of relativistic hydrodynamics are obtained in the simplest way by Landau's 
heuristic argumentation: for perfect fluid (without viscosity and heat conductivity) the energy- 
momentum tensor is 

Tf,^ ^ diag{e,p,p,p) (1) 

in the local rest frame, so we have 

Tf^ty = wui_,u^ - pgf^^ (2) 

in arbitrary frame. Projecting the energy-momentum conservation law d^T^^ = orthogo- 
nal and parallel to it^, we obtain the relativistic Euler equation and the energy conservation 
equation as 

wu^d^u^" = {g'^P - uf'uP) dpp, (3) 
wd^.u" = -ui'df.e. (4) 

The charge conservation equation (for one charge) is 

(nu^) = 0, (5) 

but if there are many different particles with the same charge, then this has to be supplemented 
by the chemical potentials. For instance, in case of baryonic or electric charge, particles and 
antiparticles carry opposite charges, and they chemical potentials are the same but of opposite 
sign, so we have 

J2 M^^M = 0. (6) 



3 Known solutions 

In this section we summarize some already known exact results: the renowned Landau-Khalat- 
nikov solution, the Hwa-Bjorken solution and other recent multi-dimensional solutions. 



3.1 The Landau-Khalatnikov solution 

Landau invented relativistic hydrodynamics in the early 1950s and proposed the application 
of it to describe high-energy particle reactions (in that time mostly cosmic ray events). In a 
collaboration with I. M. Khalatnikov, they found the first exact solution [4|5| , which is a 1-t-l 
dimensional implicit solution, with e = Sp equation of state (EoS). The Landau-Khalatnikov 
solution has nice realistic features: it is an accelerating one, the initial condition describes a 
finite piece of matter at rest, then it starts to expand. Another important feature is that the 
rapidity distribution of the particles is approximately Gaussian. However, the solution itself is 
very complicated, since the independent variables — the time and spatial coordinate — are 
given in terms of horrendous integral formulas involving Bessel functions integrated over the 
temperature and the fluid rapidity "variables" . Thus we don't quote the whole result here. 
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Table 1. The new family of solutions. 



3.2 The Hwa-Bjorken solution 

The Hwa-Bjorken solution [6|7j is a 1+1 dimensional, expanding accelerationless solution. It is 
simple and explicit, this is its main advantage compared to the Landau-Khalatnikov solution. 
It is written down in the easiest way in Rindler-coordinates r and 77, which arc defined by 

t — TCOshrj , r = Tsinh77 (7) 

with r being the spatial coordinate. The Hwa-Bjorken solution is given by the velocity and 
entropy profiles 

V = tanh rj = - , a = ag — , (8) 
t T 

where the subscript refers to the initial condition. It works for t + B = k{p + B), for arbitrary 
values of n and B. This boost-invariant solution results in a constant rapidity distribution of the 
particles. This is a rough prediction clearly not valid at the present experimental situations. 
However, Bjorken has given a simple estimate of the initial energy density of the reaction 
based on this solution, this is why it became so renowned. In the next section we present such 
solutions that are simple and explicit (as the Hwa-Bjorken solution) and are accelerating, and 
yield realistic rapidity distributions (as the Landau-Khalatnikov solution) . 



4 New simple analytic solutions 

The Landau-Khalatnikov solution contains acceleration, and yields a finite, realistic. Gaussian- 
like rapidity distribution. On the other hand, the Hwa-Bjorken solution is very simple and easy 
to handle, this is why it became so important, although the accelerationless, boost-invariant 
fiow profile and the constant rapidity distribution obviously does not agree with the results of 
real high-energy experiments. 

We found such new, analytic, explicit and simple solutions, which do not lack acceleration, 
and yield finite, realistic rapidity distributions ,l,2j . We found solutions for 1-1-1 dimensions, 
and also spherical flows in arbitrary number of spatial dimensions {d will denote their number). 
For presenting the solutions we use Rindler-coordinates, which are defined by eq. ((T]) for t > r, 
that is ,, inside the lightcone" . (From now on t is the time and r is the radial coordinate in the 
1-l-d dimensional case, and the single spatial coordinate in the 1-1-1 dimensional case.) 

For a discussion of the mathematical derivations in greater detail, apart from Refs. |l|2j . 
we recommend Ref. [3j. Here we present only the main results. The velocity and the pressure 
is given by 

v = tanHXrj) , P = Po[-) ^9) 

The value of the constants X, 4>x, d and k are constrained: different set of values yield different 
solutions. Table [1] shows the possible cases, every row of the table being a solution. In the 
following we discuss them separately. 

— Case a) isn't a new solution, it is just the well-known Hwa-Bjorken solution in l-fl di- 
mensions, and the only recently discovered Buda-Lund type of solutions |12|13j in 1-f-d 
dimensions. We quoted it for the sake of completeness, since this solution is also a member 
of the class. Apart from this case, all the other solutions possess non-vanishing acceleration. 
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— Case b) was found first, with other methods. It describes spherical expansion in d dimensions. 
The EoS is constrained in a way that k must be equal to d. (For instance, in case of three- 
dimensional expansion k — 3, which is the EoS of an ultra-relativistic ideal photon gas.) 
It is interesting to calculate the equation of trajectories R{t) from the velocity field in 
Minkowskian coordinates: 

Ofr 1 . 

"=^-72 ^ i?(t) = -(Vi + M)2 + i) (10) 

with ao — I 2^° 2 1 ■ So the fluid elements have constant uq acceleration in the local rest frame. 
ro~*o| 

— Case c) and d) were found first by T. S. Biro in 1+3 dimensions |14| . we generalized them 
to arbitrary number of spatial dimensions. These solutions are also accelerating, and since 
(/)A ^ 0, if d 7^ 1, the pressure field explicitly depends on rj, so the solution is finite in r]. 

— Case e) has a remarkably general velocity field: the A parameter can be arbitrary. We 
call A acceleration parameter, because it somehow governs the acceleration of the flow. 
On the other hand, this solution works only for d = 1 and n — 1, which is obviously 
a drawback. Nevertheless, this solution can be considered as a ,, smooth extrapolation" 
between the previous cases. In the next section we shall use this solution to calculate the 
rapidity distribution, and see that the width of it is controlled by the value of A. Hence 
in principle, the value of the parameter A of the hydrodynamical solution can be obtained 
from measurements. 

For a more thorough review of these solutions we recommend ref . [3] . 

For illustration, in Fig. [1] we plot the spatial temperature distribution at various points 
of the evolution of a collision for different A values. All the plotted solutions go to the same 
final temperature at mid-rapidity (77 = 0), although their initial temperature and acceleration 
parameter A are different. We shall see in the next section, that solutions with different accel- 
eration parameter A have different rapidity distributions, and the acceleration parameter can 
be determined from the measurement of the widths of the rapidity distribution. 



5 Rapidity distribution 

We calculated the rapidity distribution, ^ of the particles numerically and (approximately) 
analytically as well. We used the solutions discussed in the previous section as case e.), that is, 
where the A parameter is arbitrary. We assumed that at a certain freeze-out massive particles 
appear. (In the figures of this paper we used pions with m = 140 MeV, but other particles can 
be used as well ~ and their contribution can be added to the one of pions). We have chosen 
the freeze-out condition as follows: the temperature in 77 = should reach a given T/ value 
(subscript f means freeze-out) and the freeze-out hypersurface should be pseudo-orthogonal to 
the four- velocity field Uij,{x). The equation of this hypersurface is 

^ cosh ((A - 1)77) = 1. (11) 

The details of the calculation of ^ is found in refs. jl|3|. We only quote the result here: using 
a saddle-point integration in r/, for A > 1/2, m/Tf ^ 1 and Vo{s) = 1 we got 



dn dn 
Ay dy 



cosh 



(^)e-*7[™^'^°(-)-^], (12) 



where a = ■ Some typical cases are plotted in Fig. [2l This figure is just an illustration, the 
parameter values on this figure are not realistic ones. 
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Fig. 1. (Color online) These plots show the temperature distribution of diflerent one dimensional 
solutions (case e. of Table [1] with different A parameters) with different initial temperature. All the 
plotted solutions go to the same final temperature at mid-rapidity (rj = 0). The A = 1 solution 
corresponds to the Ifwa-Bjorken solution. All the A > solutions are shown within the light-cone 
(vertical solid black lines). The horizontal solid red line indicates the A = case, corresponding to a 
static, infinite, homogeneous medium with a constant temperature, and its domain extends beyond the 
liEht-cone, while for A > we plot the solutions at r < f only. 
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Fig. 2. (Color online) Normalized rapidity distributions from the new solutions in 1-fl dimensions - 
case e.) of Table[T]- for various A, Tf and m values. Thick lines show the result of numerical integration, 
thin lines the analytic approximation from eq. (|12[) . For A > 1 and not too big Tf it can be used within 
about 10 % error. 



6 Energy density estimation 

Now as we have an analytic approximation for the rapidity distribution, we are able to fit 
it to real experimental data and extract the A acceleration parameter from them. Since our 
treatment includes acceleration effects, the new solutions provide a more realistic picture of 
the initial period of high energy reactions. In this section we show how this can be used for 
improving the famous energy density estimation made by Bjorken. It is clear that initial energy 
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Fig. 3. (Color online) This figure shows that if there is no acceleration, 770 = »7/, but for the accelerating 
case, a correction factor has to be applied. 



density is a quantity of crucial importance when one wants to interpret the conclusions drawn 
from high energy experiments. 

We follow Bjorken's method and modify it when acceleration effects become important. 
Let us focus on a thin transverse piece of the produced matter at mid-rapidity, as seen on by Fig. 
2 of ref. 7 . The radius R of this slab is estimated by the radius of the colliding hadrons or nuclei: 
R = l.lS^^/'^fm. The volume is dV = (i?^7r)rd77, where r is the proper time of observation 
and dry is the space-time rapidity extent of the slab. The energy content dE = {mt)dn, where 
(mt) is the average transverse mass at mid-rapidity, so similarly to Bjorken, the initial energy 
density is 

(jntjdn 
[R-^-KjTo drio 

Here tq is the proper-time of thermalization, estimated by Bjorken as Tq w Ifm. For acceler- 
ationless, boost-invariant Hwa-Bjorken flows Vo — Vf — 2/; however, for accelerating solutions 
one has to apply a correction of 

|^|^ = (2A-l)(V.or^ (14) 

(see Fig. [3]). These two factors contain the acceleration effects on the energy density estimation, 
see ref. [3] for details. 

So the initial energy density £9 can be accessed by an advanced estimation Ec as 

£Bj \To J (R'^'kJto dy 

Here Ebj is the Bjorken estimation, which is recovered if ^ is flat (i.e. A = 1), but if A > 1, 
i.e. the flow is accelerating, both correction factors are greater than 1, so Sq is under- estimated 
by the Bjorken formula. Fig. 0] shows our fits to BRAHMS dn/dy data [T5]. From these fits we 
have found A = 1.18 ±0.01. 

Using the Bjorken estimate of esj = 5 GeV/fm'^ as given in ref. [16], and t//to = 8 ± 2 
fm/c, we find an initial energy density of Sc = (2.0 ± Q.l)eBj = 10.0 ± 0.5 GeV/fm^. 
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Fig. 4. (Color online) Left panel: dn/dy data of negative pions, as measured by the BRAHMS col- 
laboration [15] in central (0-5%) Au-|-Au collisions at ^/snn = 200 GeV, fitted with eq. (|12|l (1-1-3 
dimensional case). The fit range was —3 < j/ < 3, to exclude target and projectile rapidity region, CL 
= 0.6 %. Right panel: ec/esj ratio as a function of r//ro. 



7 Life-time estimation 

As another application, let us focus on the life-time of a high energy reaction. For a Hwa-Bjorken 
type of accelerationless, coasting longitudinal flow, Sinyukov and Makhlin [17j determined the 
longitudinal length of homogeneity as 

TBj. (16) 

Here mt is the transverse mass and tbj is the (Bjorken) freeze-out time. This result provides a 
means to determine the life-time of the reaction, if one simply identifies it with Tsj. However, 
if the flow is accelerating, the case is a little bit more complicated: the estimated origin of the 
trajectories is shifted back in proper-time, so the life-time of the reaction is under-estimated 
by TBj. From our solutions we have (for a broad but finite rapidity distribution, so that the 
saddle-point approximation can be used): 

^;ong = Y — y ^ Tc = XTBj- (17) 

Thus the new estimation of the life-time, Tc contains a A multiplication factor. BRAHMS 
rapidity distributions in Fig. |4] yield A = 1.18 ± 0.01, so they imply a 18 ±1 % increase in the 
estimated life-time of the reaction. 




8 Summary 

We presented a class of simple solutions of relativistic perfect fluid hydrodynamics. Because 
the new solutions are accelerating, but still explicit and exact, they can be utilized in the de- 
scription of high energy reactions in a more reliable way than either the Hwa-Bjorken of the 
Landau-Khalatnikov solution. We have shown, that even for a fixed equation of state, differ- 
ent initial conditions could lead to the same freeze-out distributions at mid-rapidity. However, 
these solutions of different acceleration parameter A have different width of their rapidity distri- 
butions, and the acceleration parameter can thus be determined from the measurement of the 
widths of the rapidity distribution. We have calculated the rapidity distribution of the produced 
particles, and showed two possible applications: and advanced estimation of the initial energy 
density of the high energy reaction and an advanced life-time estimation. We fitted the analytic 
approximation of the rapidity distribution to experimental data (measured by the BRAHMS 
collaboration at RHIC) and determined the A acceleration parameter of the flow. From it we 



8 



Will be inserted by the editor 



concluded that the energy density estimations based on Bjorken's method has to be corrected 
by a factor of 2. We also obtained a correction to the life-time measurements: a ~ 20% increase 
is due to the presence of acceleration in the initial evolution of the matter. Thus the effects 
of work, done even by the fluid cell at mid-rapidity, can be determined from the measurement 
of the widht of the rapidity distribution and this information is important for an advanced 
estimation of initial energy densities and the life-time of the rapidity distribution. 
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